Anisotropic conductivity of disordered 2DEGs due to spin-orbit interactions 
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We show that the disorder-averaged conductivity tensor of a disordered two-dimensional electron 
gas becomes anisotropic in the presence of both Rashba and Dresselhaus spin-orbit interactions 
(SOI). This anisotropy is a mesoscopic effect and vanishes with vanishing charge dephasing time. 
Using a diagrammatic approach including zero, one, and two-loop diagrams, we show that a con- 
sistent calculation needs to go beyond a Boltzmann equation approach. In the absence of charge 
dephasing and for zero frequency, a finite anisotropy axy oc e^ /pplh arises even for infinitesimal SOI, 
where pp is the Fermi momentum and / is the mean free path of an electron. 

PACS numbers: 72.10.-d, 72.20.-i, 71.55.Jv, 73.23.-b, 73.20.Fz, 72.25.Dc 
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I. INTRODUCTION 

The interplay between spin and charge coherence ef- 
fects in mesoscopic semiconductors produces interesting 
transport phcnomenaF^i^. They arc based on the spin- 
orbit interaction (SOI), such as Rashbaii. or Dresselhaus^ 
type, which establishes a coupling between orbital and 
spin degrees of freedom of the electron. 

One major effect of SOI on the conductivity of a dis- 
ordered semiconductor is the sign reversal of weak local- 
ization effects. In this respect the influence of the SOI is 
similar to that of a magnetic field: it increases the con- 
ductivity. For Rashba and Dresselhaus SOI and in the 
presence of a magnetic field, this "weak antilocalization" 
effect has been studied in a number of papers^i^i^. 

In the absence of magnetic fields, the Rashba SOI can- 
not violate isotropy of the energy spectrum. In this case, 
the conductivity tensor ctq,^ is invariant under rotation of 
the coordinate system (CS) by tt/2 and simultaneous sign 
change of the SOI. However, this sign is irrelevant for the 
conductivity (see Appendices), and thus axx = (^yy and 
<^xy ~ ^o'yx- In addition, due to time reversal invariance, 
o'xy = o'yx, so that aaf3 remains isotropic, i.e. a^y = 0. 
This reasoning is no longer valid if rotational invariance 
of the spectrum is broken, e.g., by a Zeeman term which 
violates both time- reversal and rotational symmetries, 
leading to a finite anisotropy of the conductivity^^. 

A similar situation arises even without magnetic 
fields but when different types of SOI (such as Rashba 
and Dresselhaus) are present: the spectrum becomes 
anisotropioii so that one can expect anisotropy of the 
conductivity. However, differently from before^^, a sys- 
tem with only SOI remains invariant under time rever- 
sal. Nevertheless, wc demonstrate below for a disordered 
two-dimensional electron gas (2DEG) that breaking of 
the rotation symmetry of the spectrum alone leads to an 
anisotropic conductivity. 

Previous calculations in such a system were based on 
the Boltzmann equatio n^^i^^ , with the outcome^ that 
the conductivity docs get enhanced by the SOI, but re- 
mains isotropic even when both Rashba and Dresselhaus 
terms are present. 

However, in the presence of phase coherence, both for 



charge and spin, the standard Boltzmann approach is 
no longer sufficient. Indeed, this approach correctly de- 
scribes contributions to aa/3 of the order of the Drude 
conductivity ctd = e'^p^l/2h, but, as is well-known, it al- 
ready 'lacks accuracy' to describe weak localization cor- 
rections ~ aY)h/p-pl <^ (7d (here pp and / are Fermi 
momentum and mean free path of the electrons, resp.). 
We will see that the required accuracy for finding the 
anisotropy of the conductivity is even higher. 

In the diagrammatic approacbi^ used below, the lead- 
ing contribution of a diagram to is of the order of 
(7d/(pfV^)"; where n is the number of loops built by 
cooperon and diffuson lines in the diagram. (This is often 
called "loop expansion" .) Thus, the zeroth order is rep- 
resented by two diagrams: the Drude "bubble" and the 
"vertex correction" , which are often referred to as "zero 
loop approximation" (ZLA)^^, indicating that these two 
diagrams have no loops made of cooperon and/or diffu- 
son lines. One can easily check that, in leading order, 
these ZLA-diagrams are SOI-independent (the "bubble" 
gives (Td, while the "vertex correction" vanishes). 

The SOI-depcndent contribution to (Jap, coming from 
the ZLA, is isotropic and on the order a-o / {p-pl / hY . We 
will see that this is of the same order as contributions 
from diagrams with two loops. Thus, a consistent calcu- 
lation (i.e. a systematic expansion in powers of h/ppl) re- 
quires consideration of all the diagrams having zero, one, 
and two loops. Below we demonstrate that the isotropic 
conductivity, obtained from the Boltzmann equatiorti^, 
corresponds to the ZLA. The inconsistency of the ZLA, 
and thus of the Boltzmann equation, has been pointed 
out beforoi^ in the context of the spin-Hall effect. 

Taking all relevant diagrams systematically into ac- 
count, wc find that the conductivity has finite anisotropic 
components given by Eq. (|16p . or, when expressed in 
terms of charge and spin dephasing times, by Eq. (|18p . 
In the fully phase coherent limit, a finite anisotropy 
axy oc /2t:pyI exists even for infinitesimal SOI. 

The paper is organized as follows. In Sec. [Hi wc cal- 
culate matrix elements of the diffuson at zero frequency; 
the result leads to the cancellation of the anomalous part 
of the velocity operator. Then, in Sec. IIIII we realize 
that the (most commonly used) zero loop approximation 
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(ZLA) gives an isotropic contribution to the conductivity 
tensor, similar to the result of the Boltzmann equationi^. 
Then, we demonstrate that this ZLA result is incomplete 
in a given order H/ppl, and higher order diagrams have 
to be considered. In Sec. IIVI we obtain the general form 
of the expansion for the anisotropic part of the conduc- 
tivity. The main results arc obtained in Sec. |Vl where 
we calculate the leading contribution to the anisotropy 
determined by two-loop diagrams and give estimates for 
real samples. 



II. KUBO FORMULA AND VERTEX 
RENORMALIZATION 

Since the considered system is invariant under time re- 
versal, (Tq/j is symmetric and can thus be diagonalized. 
We will proceed with calculations in the CS rotated by 
7r/4 with respect to the original one, where the 2D con- 
ductivity tensor is diagonal. 

In linear-response theory, the conductivity tensor is 
given by the Kubo-Grccnwood formula: 



loops composed by cooperon and diffuson lines. Each 
of these diagrams contains two velocity vertices, botbi^ 
being renormalized by the vertex correctioni^, so that 
their "anomalous" (i.e., SOI-dependent) part cancels: 



-(a + 6)(T2 
(a — 6)(Ti 
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are the components of the diffuson at zero momentum, 
K = + {xa + Xby][l + {xa — Xb)'^], and we have in- 
troduced dimensionless Rashba and Dressclhaus ampli- 
tudes, Xa = 2p-po.T/h and Xb — 2ppbT/h. [Note that Q 
is not exact at finite frequency.] 



01,(3 = x,y, 



(1) 



where the ovcrbar indicates averaging over the different 
disorder realizations, Va = j-[H, r^,] is a component of the 

velocity operator, and Gr/a ~ [Ef — H ± iO]^^ with Ep 
being the Fermi energy (the derivation of ([T]) is analogous 
to the one ioi^). The Hamiltonian in our (rotated) CS 
reads 



p 



H =^ + Vs + C/(r), 
Zra 

Vs =(a - b)cFiPy - (a + b)a2Px, 



(2) 



where a and b are the amplitudes of Rashba and Dresscl- 
haus SOI, (71^2 are expressed in terms of Pauli matrices as 
CTi,2 = (o'2 ± (Ti)/-\/2, and U{r) is a short-range impurity 
potential, U{r)U{r') = H'^{mT)~^S{r — r') (we use nota- 
tions similar to^). The conductivity ^ is not affected 
by a simultaneous sign-reversal of a and 6; thus, without 
loss of generality, we can assume that a -I- & > 0. 

In the absence of SOI, the disorder-averaged retarded 
Green function (GF) is given byi^ 
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where (Tq is the 2x2 unity matrix (due to spin degrees of 
freedom), and r is the mean time between collisions of an 
electron off impurities. The presence of the SOI modifies 
^ as follows: 



Gr = G 



(0) 



n>0 



(o)V' _ 
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The disorder-averaging of ([T]) produces an infinite num- 
ber of diagrams classified according to the number of 



III. THE ZERO-LOOP APPROXIMATION 

In the ZLA, the conductivity is given by two dia- 
grams - the "bubble" and the "vertex correction"—, 
which, when being summed, result in the following SOI- 
dcpcndent correction to the conductivityi^: 



ScTap = 17 Tr [vaGn{p)v[3GAip) 
ti p 



(0), 



xl 



h 2pFl/h VO 1 
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(7) 



where G^ is given in ([3]) and Trp = J 
suit ([7]) has been confirmed (up to a missing factor 4) 
ini^ by using the Boltzmann equation approach; so one 
can conclude that ZLA corresponds to the result of the 
Boltzmann equation with the simplest (and most com- 
mon) form of the collision integral. However, since ([7]) is 
of the order of (Jpi / [pvl / ^ ctd, the calculation in the 
ZLA (as well as the Boltzmann equationi^) is incomplete. 

To re-enforce this point, let us consider the special case 
when a = ±6. Then, all four operators under the Tr in ([T]) 
can be diagonalized by a momentum-independent unitary 
transformation, and one can sec that the SOI-dcpendcnt 
correction to the conductivity vanishes already before the 
disorder averaging (sec Appendix IB|) . This is just further 
evidence that the expression ([7]) is incomplete; other con- 
tributions must be taken into account to cancel it when 
a = ±b. 



The rc- 



IV. CONTRIBUTION OF HIGHER ORDERS 

According to the loop expansioi*i^, diagrams having 
one (weak localization^) and two loops may produce 
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FIG. 1: Three relevant two-loop diagrams (out of nine) for 
the conductivity tensor (Jq^. The matrix vertices a and a 
are written assuming clockwise (in lower Hikami boxes) and 
counterclockwise (in upper Hikami boxes) direction of writing 
traces (Tr). Each diagram contains small-momentum singu- 
larities, which cancel each other in accordance with the theo- 
rem fromi^. 



contributions to ctq/j of the same or even larger order 
than ([7]). There exist one diagram with one loop and 
nine diagrams with two loops. 

From now on we assume that the SOI and the spectrum 
anisotropy due to it are small: 



X = ^ xl + xl <C 1, (5 



2ab 



\5\ « I. (8) 



We have chosen x and 6 as expansion parameters be- 
cause they provide a uniform expansion for the conduc- 
tivity dH). 

The anisotropic part of the conductivity tensor is given 
by an expansion in inverse powers oi p-pl/h^ 1: 



'yy 



1 



(PFi/hy 



m,n>0 



, (9) 



where we used (|C6|) . 

By calculating the conductivity in the ZLA we have 
checked that = Vm, n. From the properties of the 
weak localization (WL) diagram (which is the only one 
that could contribute to S}^^), one can easily show that 



'-'oo 



as well. Thus, in the hmit when pplS'^/h < 1, 
the leading anisotropic contribution is given by the term 
oc Sqq in the expansion ([9|), which we calculate below. 



V. DIAGRAMS WITH TWO LOOPS 

We have checked that, like the ZLA, the WL diagram 
does not contribute to Sqq. The same is true for 6 (out 
of 9) two-loop diagrams. The remaining 3 diagrams, 
which we calculate below, are depicted in Fig. [1] In 
addition to the figure caption, some comments are nec- 
essary. Straight bold lines represent averaged retarded 
and advanced GFs, given by ^ and its Hermitian con- 
jugate. The squares built of GFs are Hikami boxes^^. A 



Hikami box (HB) is given by three contributions, cfJ^: 
one empty box and two boxes with an impurity line in- 
side. Wavy lines with (without) arrows represent cooper- 
ons (diffusons). 



Using the identity J2a=o^a'''^, 



= 2Ss-^s2Ss3Si, we 
write the diffuson in the basis of matrices a as Dq = 

2^ 1 — Xoiq.) , where 1 is a 4 x 4 unity matrix, 
and 



^Tr[.„G-(p).,G: 



^(p-q)]- (10) 



For cooperon lines in diagrams [TId and[Tt, we use the ex- 
pression Cq = 2^ ^ -^c(q)] , so that the cooperon 
series starts with one disorder line. Thus each diagram 
in Fig. [T] contains an extra contribution, belonging to the 
WL diagram. However, this contribution can affect 5'^i„ 
in © only for to > 1, so we neglect it. To make the 
similarity between diffuson and cooperon more evident, 
we use a different matrix basis for the cooperon: 



X, 



a/3. 



^(q-p)n, (11) 



(q) = 

where Ua 

2(5sjS2'^s3S4- Because of the time-reversal symmetry. 



CT2CTa, and ELo^'^^H 



^2GS( 



AV 

Cq = -Dq 



0^2 = Ga(p), so that Ac(q) = Xd{(\) and 
Thus the expressions for coopcron/diffuson 
lines are the same for all diagrams in Fig. [1] and the 
same is true for the lower HBs Lap-y. As a consequence, 
the sum of the three diagrams is equal to the diagram in 
Fig. [T^ with a "renormalized" upper HB Ra'p'-y' (given 
by the sum of upper HBs of all three diagrams), see the 
expression for Scr^^^ in p^ . 

Due to (jC6|) the anisotropy of the conductivity is given 
by (Jxx{a,b) - ayy{a,b) = axx{a,b) - a^xia^-b), so it 
suffices to calculate only the SOI-dependent part of axx- 
The total contribution of the diagrams in Fig. [1] is an 
integral from the product of two HBs {L and R) and 
three diffusons {D): 

h J (27r)2 J (27r)2 



a',/3',7'=0 



(12) 



First, we evaluate (|12p for zero frequency w = 0. In the 
absence of SOI, the HBs are given by 

Llp-y = -2i^^^^{kx + qx) Tr[CTc<5-/3CT7] 



lb 



Ral3f = -4z {(1 - Sf3o){kxSaoSpj + qxSyoSap) + 

[kx + 9x)(l - <5qo) [SpoSa-i + «eQ/37(l - '5/3o)(l - ^-yo)] }, 

(13) 

where Cap-y is the Levi-Civita tensor. The SOI-dependent 
part of the HBs is shown in Tab.[Tl Due to SOI, additional 
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TABLE I: The first order SOI-correction to the Hikami boxes 
(L^^^ and W^ipf^i). Other matrix elements vanish. Each of 
three upper Hikami boxes has non-zero elements with a = 
or 7 = 0, so that each diagram diverges at k ^ or at 
q — » 0; however, these elements cancel each other (note the 
zeroes in the last row), so that the complete expression for 
the conductivity converges. 



nonzero matrix elements appear, u^p^ 



that Laf}^ - L^ip^T-u^p^ 
see Tab. U 

Like the HBs, Xu defined in PT7)1 can be split into a 
main part and a SOI-dependent correction, which (for 
w = 0) are given by X^(q) = (l - Pq^/2h^) 1 and 



and Ra'/3'-y' 



and Rl,p,^,, so 
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9. 



(14) 



\0 i^{Xa+Xb) i^{Xa-Xb) xl + : 

Hence the diffusons in (fT2|) arc given by 



1-X°(q)-X|,(q; 



(15) 



where Xfy is written in (HH). 

Two angular integrations in p2|l were performed an- 
alytically using maximEfc2£ (the result is too lengthy to 
be presented here), while the integrations over fc, q were 
done numerically, using quadpack^^. 

Finally, we obtain the following result for the 
anisotropy of the conductivity at w = 0, 



r r,r II r,c2 '^XaXb e 

da = 2da^^ = a^x - Oyy = Ib^^-^ 

X., -\ 



2ttpfV (16) 



-■00 



-5.6 X 10" 



2XaXb X^ + Xb I- 



The result ((161) been obtained in the CS rotated 
by 7r/4, where the sum of Rashba and Dresselhaus terms 
are given by Vs in In the original coordinates, the 
conductivity has equal diagonal elements, and equal non- 
zero off-diagonal elements a^y = o'yx ~ 6a/2. 



The correction ()16p depends non-analytically on the 



SOI in the vicinity of 



0, so that at 







even an infinitesimally small SOI can produce a finite 
anisotropy. This is not entirely surprising, since a sim- 
ilar non-analyticity is also seen to emerge in the weak- 
localization correction with Rashba SOI only and for van- 
ishing dephasing (see, e.g., Eq. (17) in Rcfi^). 



The analyticity is restored at finite frequency: 



6a 



-2- 0.25 ■ 



—2iujT ■ 2xaXb e 1 



{xl + xf - 2iujt)^ 2Tr pfV (17) 

2xaXb x^ + xl ujT 1. 

To account for charge dephasing effects, wc follow stan- 
dard procedure^^ and substitute —iur t/t^, with 
the charge dephasing time (e.g., due to electron-electron 
interactions). In addition, one can express the SOI pa- 
rameters in terms of spin dephasing times of Dyakonov- 
Perel type^^, allowing us to rewrite p^ . (fT7)) as 

o2 



5.6 X 10 



-3 T- - T+ e 



1 



6a = 



0.13- ( 



T-f 27r Eyt ' 

2 1 



T_ I 2tT EpT ' 



T± < T0, 



T4, < T±, 



(18) 



where2^2r/r± = {xaTXh)"^- In a diffusive semiconductor 
(e.g., GaAs) for ppl/h = 5 we estimate 6a < 5 ■ lO^^tro 
in case t± T,f, and 6a < IO^^ctd hi the opposite case 
T± > r^. 

Thus, in a system with finite dephasing time T0, the 
anisotropic conductivity (|18p is analytic in t±. However, 
in a fully phase-coherent system with = 00 (which is 
commonly predicted at zero temperature), the conduc- 
tivity tensor depends non-analytically on the SOI. 



VI. CONCLUSION 

The combination of Rashba and Dresselhaus SOI leads 
to anisotropy in the energy spectrum, which, in turn, 
leads to anisotropy in the conductivity of a disordered 
2DEG. This anisotropy, being due to phase coherence of 
charge and spin, comes from two-loop diagrams in the 
perturbative approach, while it is absent in the Boltz- 
mann equation approach. In the limit of full phase co- 
herence, the system "jumps" from isotropic to anisotropic 
state for infinitesimally small SOIs. 

Wc arc grateful to B. Altshulcr, M. Duckheim, 
A. Khaetskii, and E. Rashba for helpful discussions. We 
acknowledge financial support from the Swiss NSF and 
the NCCR Nanoscience. 



APPENDIX A: CONDUCTIVITY FOR 
PURE-RASHBA OR PURE-DRESSELHAUS SOI 

Using the time-reversal invariance we argued in the 
Introduction (Sec. [l| that the conductivity tensor is 
isotropic in case of a purc-Rashba SOI. Here we demon- 
strate this isotropy for both pure-Rashba and pure- 
Dresselhaus cases using only rotation symmetry of the 
disorder-free part of the Hamiltonian. 

The conductivity tensor in the rotated CS is connected 
with its value in the original CS through the transforma- 
tion 



0^ 



(Al) 
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where is a SD-rotation matrix in the (x, ?/)-plane by 
angle (j). In this section, we prove that when a = or 
6 = 0; the conductivity tensor is invariant with respect 
to rotations in the (x, ?/)-plane; together with the require- 
ment (|Aip this means that aafs is isotropic. 
We define matrices 



(A2) 





The Hamiltonian ^ can be rewritten as 

^2 



H 



P 
2m 



+ aVR + bVD + Uir), 



where (in the original -unrotated- CS) 

Vr = e, • [cr X p] , Vd = • [Ccr x p] 



(A3) 



(A4) 



cr = (cTi, (T2, (73). and is the unit vector along the z- 
axis. A rotation of the CS corresponds to the transfor- 
mation of coordinates, momenta and spins according to 

(fJo, (y) {aa,R4,cr), p R^p, r R^r. (A5) 



from any other diagram is invariant under the transfor- 
mation (|A5p in case when 6 = 0, so that [due to lAl[) ] the 
pure-Rashba conductivity tensor is isotropic. 

Let us now demonstrate that also in case of pure- 
Dresselhaus SOI the conductivity is isotropic. The Dres- 
sclhaus SOI can be transformed into Rashba SOI by 
a unitary transformation: Vd = T^VrT; the same is 
true for the pure-Dresselhaus velocity operator: v^^ = 

#^+5Vn,r„ = T'^v^^T. So, one can transform 

the pure-Dresselhaus conductivity tensor into the pure- 
Rashba one: 



G 



Dbx-.DbQDbl 



R ''^IS '-'A J 



(A8) 



so that af^ 



(J^p. With (|A8p we proved this statement 



for the (simplest) bubble diagrami^; it can be general- 
ized to higher-loop corrections analogously to the pure- 
Rashba case [sec (|A7|) above]. We conclude that the 
disorder- averaged conductivity tensor is isotropic in a 
system with either pure Rashba or pure Dresselhaus SOI. 



We note that Vr is invariant under the transformation 
[1: 



(A6) 



Vr = • [i?0(T X i?0p] = 




The same invariance holds for the (disorder-averaged) 
pure-Rashba Green functions G^"^. Thus, the contribu- 
tion to the conductivity tensor from the (simplest) bubble 
diagranii^ 



= Tr Tr [v 

fl p spin 
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2m 



R ''[3 

aVR,ra 



G \ 



(A7) 



is invariant under the rotation (jASp . and hence isotropic 
[due to lAip ] in a pure-Rashba SOI system. As an exam- 
ple of a higher-loop correction, let us consider the two- 
loop contribution (jl2p . The symmetry transformation 
(|A5p will affect ([T^ twofold: (i) as a momentum rotation 
k, q — !■ i?0k, i?0q and (ii) through a different set of Pauh 
matrices ((Tq, cr') ~ (do, RtfiCr) in the calculation of diffu- 
sons and Hikami boxes. The rotation of momenta can be 
absorbed by a variable change in the integration. Then, 
the result of the calculation does not depend on the choice 
of Pauli matrices, once they obey standard commutation 
relations. Since the new set (ctq, cr') is connected through 
a unitary transformation with the standard one, (cto,"") 
the commutation relations are preserved so that (|12p is 
invariant under the transformation (|A5p . Thus, the two- 
loop contribution IT^ does not change when the CS is 
rotated in the [x, y)-plane. Analogously, the contribution 



APPENDIX B: CONDUCTIVITY IN CASE 
WHEN a^±h 

The case, when the moduli of Rashba and Dresselhaus 
SOI amplitudes are equal, is speciali^i^i^. Both Hamil- 
tonian and velocity operators can then be diagonalized 
in spin space by unitary transformations: for a = 6, up 
to an irrelevant constant (in the rotated CS) 



UiHUl 



Px_ 
771 



2aa:: 



C/(r), (Bl) 



and 



where the transformation matrix is given by 
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777 / 777^ 



and 



C/2VC/J 
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+ U{r), (B4) 
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APPENDIX C: FURTHER SYMMETRIES OF cr^^ 



U2 = [uir' = 



1 —er^^/^ 



\/2 ll 



(B6) 



Let us denote the (unaveraged) Green functions without 
SOI {a = b = 0) as g,/ 

With [/i 2 from (jB3p and (jB6p we diagonahze aU oper- 
ators under the trace Tr in the conductivity expressions 
below. Then the spin part of Tr = Trgpin Trp sphts into 
two terms. For a = h, 



While the Rashba SOI Vr is invariant under rotations 
[sec (|A6|) ]. the Dresselhaus SOI Vd docs not possess this 
symmetry - rotation by an angle (j) transforms it as fol- 
lows: 



a = 6 ^ Tr Tr 

p spin 



= Tr Tr 

p spin 
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Vd ^ e^- [CR^a x iJ^p] , 
which is in general different from 
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= Tr 
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= 2Tr 
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Px . Px . 
5r ffa 

m m 
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TO TO 
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—gi—9a 

TO TO 



In the special case of rotation by = 7r/2, the anticom- 
mutator 
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PynPyn 

TO TO 



(B7) 



-b Tr Tr 

p spin 



and for a = —b, 

VyGni'yGA 
U2VyUlU2GR ul U2 vy ul C/2 G A ul 



/O \ 

{i?,/2,G} = (C3) 
\0 -2/ 

is zero in the (x, ?;)-plane, so that the r.h.s. of (jCip differs 
from (jC2p only by the sign. On the other hand, 
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Thus, we see that in the rotated CS Gaa are the same as 
in the absence of SOI for a = ±&: 



= ±b 



Uxx = o-yy = <Jxx{a = b = 0). (B9) 



so that rotation by tt/2 is equivalent to the sign change 
of the Dresselhaus SOI amplitude b, and 
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